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ABSTRACT 

We have analyzed XMM-Newton Optical Monitor (OM) UV (180-400 nm) 
data for a sample of 33 galaxies. 30 are cluster member galaxies, and nine of 
these are central cluster galaxies (CCGs) in cooling flow clusters having mass 
deposition rates which span a range of 8 — 525 Mq yr~^. By comparing the ratio 
of UV to 2MASS J band fluxes, we find a significant UV excess in many, but 
not all, cooling flow CCGs, a finding consistent with the outcome of previous 
studies based on optical imaging data (McNamara & O'Connell 1989; Cardiel, 
Gorgas, & Aragon-Salamanca 1998; Crawford et al. 1999). This UV excess is a 
direct indication of the presence of young massive stars, and therefore recent star 
formation, in these galaxies. Using the Starburst99 spectral energy distribution 
(SED) model of continuous star formation over a 900 Myr period, we derive star 
formation rates of 0.2 — 219 Mq yr~^ for the cooling flow sample. For 2/3 of this 
sample it is possible to equate Chandra/XMM cooling flow mass deposition rates 
with UV inferred star formation rates, for a combination of starburst lifetime and 
IMF slope. This is a pilot study of the well populated XMM UV cluster archive 
and a more extensive follow up study is currently underway. 

Subject headings: cooling flows — galaxies:clusters:general — galaxies:elliptical and 
lenticular, cD — galaxies: stellar content — stars:formation — ultraviolet: stars 



1. Introduction 

The 'mystery' of cooling flows (Fabian et al. 1984; Sarazin 1986; Fabian 1994) has baffled 
astronomers for over 20 years. Though the X-ray temperature often drops by a factor of three 
within the central 100 kpc (Allen et al. 2001), there is seldom enough mass at temperatures 
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cooler than ~ 2 keV to be consistent with constant pressure coohng flow models. Recent 
XMM and Chandra results show that the distribution of temperature vs. emission measure 
also does not match the prediction of these models (Peterson et al. 2003; Voigt & Fabian 
2004). While the XMM and Chandra results have reduced the implied cooling rate in most 
clusters there is still an apparent discrepancy of an order or magnitude or more between the 
X-ray cooling rates and the observed star formation rates (SFR) in many clusters. Contrary 
to conventional wisdom that there is no star formation at all in the central cluster galaxies, 
a number of contemporary studies have found significant star formation (Melnick, Gopal- 
Krishna, & Terlevich 1997; Smith et al. 1997; Hansen, Jorgensen, & Norgaard- Nielsen 1998; 
Mittaz et al. 2001; McNamara et al. 2004a; O'Dea et al. 2004; McNamara 2004b). However 
we are still far from a complete reconciliation between X-ray mass deposition rates and 
alternatively derived star formation rates in the majority of cooling flow clusters. 

While many cluster cooling flows have been well observed in the radio and with emission 
line studies (Heckman et al. 1989, and others) there have been very few published studies 
of the UV properties of cooling flow clusters. In this study we use the readily available 
archival XMM Optical Monitor UV observations to derive the UV luminosities of a moderate 
sample of clusters. Since the dominant contribution to UV flux comes from short-lived main 
sequence stars, and the total number of these stars is directly proportional to the star 
formation rate (Kennicutt 1998) the UV is one of the prime routes to understanding star 
formation. We attempt to address and quantify the connection between cooling flows and 
star formation in our pilot sample. Unless otherwise noted, this paper assumes a cosmology 
of Ho = 70 km s"^ Mpc'S Qm = 0.3, and Qa = 0.7 



2. The XMM-Newton Optical Monitor 

The XMM Optical/UV Monitor Telescope is a modified 30cm Ritchey-Chretien tele- 
scope mounted alongside the X-ray mirror modules (Mason et al. 2001) which allows simul- 
taneous X-ray and UV/optical observations of XMM targets within the central 17 square 
arcminute region of the X-ray field of view. The 256 x 256 pixels have an angular resolution 

of ~ 1" . The sensitivity limit of the Optical Monitor (OM) is 24 mag. This study primarily 
utilizes the UVWl fiher (220-400 nm), and the UVW2 filter (180-260 nm), and the UVM2 
filter (200-280 nm) for calibration purposes. Previous XMM OM studies of clusters have 
been limited to the detection of UV mission in Abell 1795 (Mittaz et al. 2001). 
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3. Sample and Observations 



Our pilot sample consists of 33 galaxies and eleven clusters. Three of them (M49, M60, 
and M86) were selected because they possess lUE photometric information (Burstein et al. 
1988) and were used to compare OM to lUE fluxes. 21 galaxies in the IC 1860 and Coma 
fields were used to calibrate the relationship between UV and IR luminosities in passively 
evolving early type and cD galaxies. Two UVWl observations of Coma were included, one 
centered on the cD galaxies, and one which is off-center and contains numerous early-type 
galaxy detections within the viewing field. 

The remaining nine clusters in our sample were chosen because of significant X-ray de- 
termined mass deposition rates, general overlap with the samples of McNamara & O'Connell 
(1989) and Cardiel, Gorgas, & Aragon-Salamanca (1998), moderate to low redshift, inclu- 
sion in the 2MASS galaxy sample, and having archival XMM UVWl data. The majority of 
clusters in our sample are low redshift objects {z < 0.1), but we have also included three 
moderate redshift clusters {z ~ 0.25) which are known to possess very strong cooling flows. 

Figure 1 shows one of the coohng flow clusters in our sample, with XMM-Newton UVWl 
contours overlaid. Table 1 lists the objects in our current sample, their redshifts, relevant 
filters, and OM exposure times. 



We fitted the flux in concentric circles around the emission peak of each galaxy to a 
radial surface brightness profile of the form: 



Where Iq is the normalization, tq is the core radius, and Ib is a constant representing the 
intensity contribution of the background. 

The surface brightness profiles are well fit by these models, cf. Figure 2. The fit 
parameters are in Table 2. 



4. Surface Brightness 




(1) 
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5. Photometry 

We retrieved pipeline-processed UV data from the XMM-Newton Science Archive and 
combined observations with the same pointing. Fhix determinations for our cahbration 
sample used a circular aperture of 7" radius, a good approximation of the lUE aperture 
(Burstein et al. 1988), and the size of the fixed 2MASS photometric aperture. Photon 
counts in the aperture were corrected for background using a 35" radius background region 
taken from an outlying area of the chip in the same observation. Final fluxes were determined 
by employing the OM counts-to-flux conversions, and errors were assessed via -\/N photon 
statistics. 

The results from UVW2 and UVM2 observations of M49, M60, and M86 were then 
compared with the lUE determined flux value of each object (Burstein et al. 1988) for the 

bandpass of the appropriate filter. Our flux measurements are consistent with theirs to within 
< 35%, errors which are in keeping with those of Vassiliadis et al (1996) and Koratkar et al. 
(1997). Two of the three galaxies have lower OM fluxes than lUE fluxes, indicating that the 
Optical Monitor does not have a red leak. 

We then determined UVWl fluxes for the remainder of our sample (see Table 1), for 
a 7" radius region, in the method described above. The UVWl data were corrected for 
average galactic extinctions (Cardelli, Clayton, & Mathis 1989). J band extinction values 
were retrieved from the NASA/IPAC Extragalactic Database (NED). 

6. Calibration 

The low redshift {z ~ 0.023) calibration sample has 23 galaxies in IC 1860 and Coma 
with 2MASS J band photometric measurements. Two of these galaxies exhibited high UV 
fluxes with respect to the rest of the calibration sample, and were described as star forming 
in a FOCA study Donas, Millard, & Laget (1995) and thus are not appropriate for inclusion 
in our UV/IR luminosity calibration of passively evolving early type galaxies, and were 
removed from the remainder of this study. 

A least squares fit between the UVWl and J band luminosities for the remaining 21 
early type galaxies is well described by 



Luv — Ci X Ljj^^ (2) 
Where Luv is the UVWl luminosity and Lir is the 2MASS J band luminosity, inside 
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7" radius apertures with Ci = 0.090 ± 0.007 and C2 = 0.708 ± 0.007 respectively (Figure 
3). We use this correlation to predict the UV light from non-starforming, passive galaxies 
assuming that the UV light from such objects is produced by old evolved stars coming from 
the same population that produces the IR light. 

7. Photometric Results 

The UV luminosity excesses of the cD galaxies in our nine coohng flow clusters were 
calculated by subtracting the expected UV light due to old stars (based on the fit obtained 
in Section 6) from the measured UV luminosity of each cluster. Figure 4 shows the J band 
luminosity of each galaxy in our sample plotted against its UV/IR luminosity ratio. Wc have 
not attempted to estimate and correct for internal dust absorption, and thus these excesses 
provide a lower limit on the UV emission in these clusters. 

Because we are analyzing data from a single filter, it is difficult to break the degeneracy 
between age and mass of the young stellar populations. We folded through the UVWl filter 
Starburst99 (Lcithcrcr 1999) rcdshifted models corresponding to continuous star formation 
over a a variety of timescales with solar abundance (Z = 0.020) and IMF powerlaw indices 
of both 2.35 and 3.3 to estabhsh the relation between OM UV counting rates and star 
formation rate. The measured UV luminosity excesses of the clusters were then used to 
estimate per-year star formation rates (Table 3) that were then compared to X-ray derived 
mass deposition rates from Voigt & Fabian (2004), Peterson et al. (2003), and R. Miishotzky 
(2005, private communication). Also in Table 3 are estimates of the duration of continuous 
star formation, and best fitting IMF powerlaw index. These values were derived by plotting 
the relationship between star formation duration and rate for each given UV excess and IMF 
slope, then finding the set of parameters (slope and duration) at which X-ray determined 
mass deposition rates intercepted those curves. The 900 Myr models should be closest to 
the average hfetime of a cooling flow (Allen et al. 2001). 

8. Summary and Discussion 

We find a significant UV luminosity excess in the central galaxies of most cooling fiow 
clusters. This implies that star formation has recently been or is currently occurring most 
likely as a direct result of radiative coohng of the intracluster medium (ICM). These con- 
clusions are qualitatively consistent with the findings of previous optical imaging studies of 
star formation in cooling flow clusters (McNamara & O'Connell 1989; Cardiel, Gorgas, & 
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Aragon-Salamanca 1998; Crawford et al. 1999). 

In the absence of multiwavelength spectral information, it is difficult to constrain a 
particular model of the star formation history of these galaxies. However we note that for 
2/3 of our sample it is possible to find some combination of the parameters of IMF slope and 
starburst lifetime which allow us to equate X-ray determined mass deposition rates with UV 
inferred star formation rates. 

In none of the OM, Chandra or HST images is there evidence of a significant excess 
central flux, therefore we expect the AGN UV contribution to our signal is small. In addition, 
we have investigated the amount of UV flux which might come from the cooling flow itself, 
and have found it to be much smaller than measured UV excesses. Since the calibration 
sample consists entirely of low redshift early type galaxies, there is a possible systematic 
error in our calculations of the SFRs of the three moderate redshift {z ~ 0.25) clusters in 
our sample, due to evolution in the UV luminosity of the old stellar population (Brown et 
al. 2003) and/or uncertainties in the redshift evolution of the star forming component. In 
future work we will attempt to calibrate a higher redshift passive E galaxy sample. 

While there is significant UV flux in 2/3 of our cooling flow sample, 1/3 of the objects 
with significant cooling rates from modern data do not have detectable UV excesses, and 
the ratio of UV light to X-ray coohng rate (Figure 5) shows a wide range. The origins of 
this scatter are unclear and a larger UV sample, with direct comparison of X-ray and radio 
imaging and IR data, is necessary to get a clearer picture of the origin of the UV hght. The 
discovery of copious UV emission from several cooling flow clusters confirms the previous 
work and shows that a significant fraction of the cooling gas in clusters can indeed form 
stars. 

Our star formation estimates are very basic, being based solely on a simple conversion 
from UV light to SFR (Kennicutt 1998). We have bypassed issues of the age of the star 
formation, the effects of reddening, the slope of the IMF and metallicity effects. The fact 
that we can, roughly, obtain similar values of X-ray inferred cooling rate and UV inferred 
star formation rate indicates that, given the above degeneracies, more and better data (such 
as could be obtained by Galex or HST spectroscopy) is necessary to determine precise star 
formation rates. 

This paper is a pilot study and we are currently undertaking an expansion of this 
project toward higher redshift, utilizing more of the OM UV filters and including the many 
additional clusters that remain in the archive. Through this expansion we hope to achieve a 
more thorough characterization of recent star formation in cooling flow clusters of galaxies. 
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image, and have count values of 14, 25, 44.75, and 80. The large splotch of UV contours to 
the lower right of the cD galaxy is due to scattered hght. 
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Fig. 2. — Surface Brightness Profile. Radial surface brightness profile of Abell 1795 
made from UVWl (220-400 nm) data accumulated in ~ 1" annular bins. A solid line traces 
the best /3 model fit. The horizontal dotted line represents the best fitting background value. 
The best fitting /3 model has a core radius of 4.09 " and beta value of 0.703 



- 11 - 



J Band Luminosity vs. UV Lunninosity 
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Fig. 3. — UV vs. IR luminosity for non-cooling flow cluster galaxies in our sample. The line 
indicates the best fitting relationship between UV and IR luminosity for passively evolving 
and non-cooling flow cD galaxies, represented by squares and triangles, respectively. Dia- 
monds indicate the two star forming galaxies (Donas, Millard, & Laget 1995) which were 
excluded from the fit. A relationship between UV and IR luminosity is expected here, be- 
cause passively evolving cluster galaxies have similar spectral energy distributions (SEDs). 
The similar redshifts of IC1860 and Coma allow us to fit this relationship without making 
distance corrections to the data. 
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J Bond Luminosity vs. UV/IR Ratio 
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Fig. 4. — UV/IR luminosity ratio of selected galaxies vs. their 2MASS J band luminosities. 
All luminosities were calculated in 7" radius regions, and corrected for galactic extinction. 
Diamonds indicate cooling flow cluster cDs, all of which are at low redshift {z < 0.1) except 
for the three points farthest to the right, which have redshifts of ~ 0.25. The line designates 
the best fitting relationship between UV and IR luminosity for passively evolving and non- 
cooling flow cD galaxies (Section 6). The high UV/IR ratio of the cooling flow clusters is a 
strong indication of the existence of hot young stars and thus a direct indication of copious 
recent star formation. 
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Fig. 5. — Excess UV luminosity vs. M. Luminosity excesses were determined by subtracting 
the best fit ratio of UV liglit due to old stars (Figure 2) from the measured UV luminosity 
of each cluster. X-ray mass deposition rates were taken from recent papers (Peterson et al. 
2003; Voigt & Fabian 2004). 
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Table 1. XMM- Newton Optical Monitor Observations 



Object 


z 


Filter 


Exposure [s] 


IC 1860 


0.0229 


UVWl 


5000 


PKS 0745-191 


0.1028 


UVWl 


3999 


Hydra A 


0.0538 


UVWl 


3000 


Zw3146 


0.2906 


UVWl 


4000 


NGC 4406 (M86) 


-0.0008 


UVWl 


3600 






UVM2 


3500 


NGC 4472 (M49) 


0.0033 


UVWl 


2000 






UVM2 


1000 


NGC 4649 (M60) 


0.0037 


UVWl 


2499 






UVW2 


2900 


Coma (center) 


0.0231 


UVWl 


3900 


Coma (outskirts) 




UVWl 


1498 


Abell 1795 


0.0633 


UVWl 


2499 


Abell 1835 


0.2520 


UVWl 


2000 


MS 1455.0+2232 


0.2578 


UVWl 


2000 


Abell 2029 


0.0779 


UVWl 


1200 


Abell 2052 


0.0345 


UVWl 


2000 


MKW3s 


0.0453 


UVWl 


2002 



- 15 - 



Table 2. Beta- Model Fits 



Object 


Rcore [ ] 


P 




T a 
-'B 


Reduced 


PKS 0745-191 


2_83±o.o4 




625.3±i2.4 


481.0^1-6 


0.691 


Hydra A 


3449±o.ooi 


0_8387±o.ooo2 


4829.1^2.4 


356.59^°°^ 


4.200 


Zw3146 


g_2i±o.oi 


2,o54±o-oo7 


1119.2^^-3 


413.2^0-^ 


0.947 


Abell 1795 


4.09±o-oi 


0,703±o.oo2 


2334.8^13.4 


336.6^°-3 


0.946 


Abcll 1835 


3.65±o-4i 


^ 253±o.i65 


288o.7±i53.o 


349.6^1-s 


0.514 


MS1455.0+2232 


4.99±o.o3 




446.9±3-5 


337.8^°-^ 


0.657 


Abell 2029 


3.68±°-02 


Q 499±o.ooi 


613.2±3-5 


449.7±o.6 


1.101 


Abell 2052 


16g±0.28 


0402±o.O22 


1593. l^i^'^-^ 


426.2±'^-^ 


1.912 


MKWSs 


128±o-33 


352±o.o22 


1018.0±i38-9 


398.6±^-i 


1.110 



^Surface brightness / in units of 10 ergs sec ^ cm ^ arcsec 
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Tablc 3. 



Cluster 


z 


UV Excess 


M2.35 " 


M3.3 " 


M (X-ray) 




a 






[10^3 erg/s] 


[MeM 


[Mo/yr] 


[MeM 


[Myr] 




PKS 0745-191 


0.1028 


11.94±o-7 


129±^ 


237±i3 


317+1 ^ 


200 


3.3 


Hydra A 


0.05384 


0.88±o-02 


9.5±o-2 


17.5±o-4 


8ir 


900+ 


2.35 


Zw3146 


0.2906 


847±o.4 


91±^ 


168±8 


525lgQ 


11 


3.3 


Abell 1795 


0.06326 


0.85±°-02 


9.1±o-2 


lg7±0.4 


30t^^ 


100 


3.3 


Abell 1835 


0.251982 


11.41±o-4 


123^5 


226^9 


Q/i+43 b 


900+ 


2.35 


MS 1455.0+2232 


0.2578 


3.36±°-^ 


36±^ 


66±« 


2001?° 


150 


3.3 


Abell 2029 


0.07795 


O.i8±o.o5 


2 y±0.5 


32±o.9 


30li° b 


7 


3.3 


Abell 2052 


0.03446 


0_19±0.007 


i9±o.o8 


0.4±o.i 


lOOl^Oc 








0.01531 


0.01™ 


o.i±"-i 


0.2±"-2 


1 r + 10 c 







^Assuming constant star formation over 900 Myr. The subscript indicates the IMF power- 
law slope of the model used. 

''Voigt & Fabian (2004). 

^Peterson et al. (2003). 



